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Human coronaviruses (HCoVs), including SARS-CoV and MERS-CoV, are zoo- 
notic pathogens that originated in wild animals. HCoVs have large genomes that 
encode a fixed array of structural and nonstructural components, as well as a 
variety of accessory proteins that differ in number and sequence even among 
closely related CoVs. Thus, in addition to recombination and mutation, HCoV 
genomes evolve through gene gains and losses. In this review we summarize 
recent findings on the molecular evolution of HCoV genomes, with special 
attention to recombination and adaptive events that generated new viral species 
and contributed to host shifts and to HCoV emergence. 


Human Coronaviruses Are Zoonotic Pathogens 

The recent emergence of severe acute respiratory syndrome-related coronavirus (SARS-CoV) 
and of Middle East respiratory syndrome-related Coronavirus (MERS-CoV) (order Nidovirales, 
family Coronaviridae, subfamily Coronavirinae) as dangerous zoonoses stirred great interest in 
the ecology and evolution of coronaviruses. Before the SARS-CoV epidemic only two HCoVs 
were known: HCoV-229E and HCoV-OC48. Two additional HCoVs, HCoV-NL63 and HCoV- 
HKU1, were discovered in 2004-2005 from clinical specimens [1]. These viruses originated in 
animals and are mainly responsible for respiratory diseases in humans (Figure 1A, Key Figure). 
Specifically, all HCoVs are thought to have a bat origin, with the exception of lineage A beta- 
CoVs, which may have reservoirs in rodents [2]. The phylogenetic relationships of HCoVs and 
other animal CoVs mentioned in this review are summarized in Figure 1A. 


A number of field studies identified and sequenced viruses related to HCoVs in wildlife 
reservoirs, and phylogenetic reconstruction provided important clues on the most likely 
events that led to the introduction of HCoVs in human populations. Several recent excellent 
reviews delve into the knowns and unknowns of HCoV origin in terms of reservoir species, 
amplification host, and, more generally, of CoV ecology [1,3-5]. In this review we instead 
focus on the molecular evolution of HCoV genomes. The general concepts of evolutionary 
analyses in viruses are outlined in Box 1, whereas the most common approaches that were 
applied to the analysis of CoV sequence evolution in terms of phylogenetic reconstruction, 
detection of recombination, and identification of selection signatures are summarized in 
Boxes 1 and 2. 


HCoV Genome Organization 

CoVs are positive-sense, single-strand RNA viruses with a likely ancient origin, and HCoVs 
repeatedly emerged during the past 1000 years (Box 3). All CoVs have nonsegmented genomes 
that share a similar organization. About two thirds of the genome consists of two large over- 
lapping open reading frames (ORF 1a and ORF 1b; see Glossary), that are translated into the 
pptaand pp1ab polyproteins. These are processed to generate 16 nonstructural proteins (nsp1 
to 16). The remaining portion of the genome includes OFFs for the structural proteins: spike (S), 
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Human coronaviruses (HCoVs) are zoo- 
notic pathogens with large and complex 
genomes. Some HCoV accessory pro- 
teins were acquired from host genes, 
and some were lost or split during 
HCoV evolution. Most likely SARS- 
CoV ORF8 became dispensable during 
the shift to the human/civet host. 


HCoV spike proteins adapted to use 
diverse cellular receptors. This 
occurred by divergence followed, in 
some cases, by convergent evolution 
to bind the same receptor. 


Recombination and positive selection 
shaped the diversity of CoV genomes, 
especially the S gene. Positive selec- 
tion in the S gene of MERS-CoV and 
related CoVs mainly acted on the hep- 
tad repeats. 


In MERS-CoV and other lineage C 
beta-CoVs, positive selection targeted 
the nonstructural components, parti- 
cularly ORF1a. Most adaptive events 
occurred in nsp3, which acts as a viral 
protease and contributes to suppres- 
sion of interferon responses. 
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Box 1. Molecular Evolution in Viruses: General Concepts 


RNA viruses are rapidly evolving pathogens that can accumulate considerable genetic diversity in relatively short time 
periods. This is mostly due to their high nucleotide mutation rates (but see text for CoVs). The diversity of extant viral 
sequences can be analyzed to construct phylogenetic relationships among species/strains and to infer the underlying 
evolutionary patterns. In the presence of recombination a single phylogenetic tree is unable to describe the evolution o 
homologous sequences. Because recombination is common in many viruses, including HCoVs, the evolution of viral 
genomes is best modeled by several phylogenies, one for each nonrecombinant fragment (Box 2). By reassorting 
mutations, recombination has the potential to generate novel viral phenotypes. Thus, not only recombination is of interest 
per se, but failure to account for its presence can distort phylogeny-based analyses, including estimates of natural 
selection [66]. Natural selection acts pervasively on viral sequences. When coding regions are concerned, natural 
selection is commonly estimated in terms of w (also referred to as dN/dS) — that is, the observed number o 
nonsynonymous differences per nonsynonymous site (AN) over the observed number of synonymous differences 
per synonymous site (dS). Under neutral evolution, w is expected to be equal to 1, as the rate at which amino acid 
substitutions accumulate is similar to the rate for synonymous changes. Due to the fact that essential protein domains 
can often tolerate only minor sequence changes, most amino acid replacements are eliminated by selection; this 
generates w values <1, a situation referred to as negative (or purifying) selection. Nevertheless, amino acid replacements 
can be advantageous for a virus in terms, for example, of host adaptation or immune evasion: in this case w values can be 
higher than 1 (positive selection). Thus, evaluation of how w varies from site to site or from branch to branch in a 
phylogeny is commonly used to describe selective events. Some possible caveats should nevertheless be kept in mind. (i) 
The saturation of substitution rates (especially dS) may occur and affect evolutionary inference when fast-evolving 
sequences are analyzed (see Box 3 for an example, and Box 2 for methods to overcome this problem). (ii) In viral 
genomes synonymous substitutions are not always neutral; this may be due to the presence of overlapping reading 
frames, conserved RNA secondary structures, packaging signals, and other functional elements (Box 3). (iii) A relaxation 
in the intensity of both negative and positive selection may occasionally occur (Box 2 and Figure 2A). 


envelope (E), membrane (M) and nucleoprotein (N). A variable number of accessory proteins are 
also encoded by distinct viruses (Figure 1B). 


Among RNA viruses, CoVs have exceptionally long genomes (up to 32 kb). Genome expansion 
in CoVs is believed to be at least partially mediated by increased replication fidelity. Although 
estimates of the mutation rate for CoVs differ, possibly depending on the phase of CoV 
adaptation to novel hosts, several studies have shown that these viruses may possess an 
unusually high replication fidelity [6-8]. Indeed, a major step that allowed genome expansion in 
CoVs and, more generally, in Nidovirales, was the acquisition of a set of RNA-processing 
enzymes that improved the low fidelity of RNA replication [9]. These enzymes include an RNA 3’- 
to-5’ exoribonuclease (ExoN) and possibly an endoribonuclease (NendoU) [9]. Additional evi- 
dence, though, suggests that features distinct from replication fidelity underlie genome expan- 
sion in Nidovirales. These include a peculiar genome organization [9] and a processive replication 
complex [10]. 


Importantly, CoV genome expansion allowed the acquisition and maintenance of genes encod- 
ing diverse accessory proteins that may promote virus adaptation to specific hosts and often 
contribute to the suppression of immune responses, as well as to virulence. Accessory proteins 
differ in number and sequence even among CoVs belonging to the same lineage (Figure 1B), 
raising interesting questions about their origin and evolution. 


Gene Gains and Gene Losses 

The acquisition (or loss) of novel protein-coding genes has the potential to drastically modify viral 
phenotypes. Thus, tracing these gain/loss events may identify important turning points in viral 
evolution. 


Among SARS-CoV accessory proteins, the origin of ORF8 has remained mysterious for a while, 
as SARS-CoV-related (SARSr) bat viruses were isolated but found to encode divergent ORF8 
proteins (amino acid identity with SARS-CoV ORF8 around 33%) [1 1-138]. Very recently, SARSr- 
BatCoVs from Rhinolophus sinicus (Rs) and Rhinolophus ferrumequinum (Rf) were isolated 
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dN: the observed number of 
nonsynonymous substitutions per 
nonsynonymous site. 

dS: the observed number of 
synonymous substitutions per 
synonymous site. 
Hemagglutinin-esterases (HEs): a 
family of viral proteins that mediate 
binding to O-acetylated sialic acids. 
Homology: the relationship between 
elements (€.g., genes, proteins) 
deriving from a common ancestor. 
Lectins: a group of proteins with 
carbohydrate recognition activity. 
Lectins are categorized in many 
distinct families depending on 
structural and functional properties. 
Maximum likelihood (ML): is a 
statistical method for estimating 
population parameters from a data 
sample. Given one or more unknown 
parameters and a sample data, the 
ML estimates of the parameters are 
the values maximizing the probability 
of obtaining the observed data. 

Open reading frame (ORF): the 
part of a reading frame that contains 
no stop codons. An ORF is a 
continuous stretch of nucleotide 
triplets that have the potential to 
code for a protein or a peptide. 
Phosphodiesterases (PDEs): are 
enzymes that break a phosphodiester 
bond. PDEs belonging to the 2H 
family are characterized by two H-®- 
[S/T]-® motifs (where ® is a 
hydrophobic residue) separated by 
an average of 80 residues. 

Positive selection: the accumulation 
of favorable amino acid-replacing 
substitutions, which results in more 
nonsynonymous changes than 
expected under neutrality (AN/dS>1). 
Purifying selection: the elimination 
of deleterious amino acid-replacing 
substitutions, which results in fewer 
nonsynonymous changes than 
expected under neutrality (AN/dS < 1) 
(it is also referred to as negative 
selection). 

Viroporins: hydrophobic viral 
proteins that can promote the 
formation of channels following 
insertion into the host cell membrane 
and oligomerization. 
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Key Figure 
Phylogenetic Relationships and Genome Organization of Human and Animal Coronaviruses (CoVs) 
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Figure 1. CoVs that infect nonhuman mammals are included only if they are mentioned in the text for comparative purposes. (A) The phylogenetic tree of complete 
genome sequences of HCoVs and selected mammalian CoVs was obtained with RAXML 8.2.4 [68]. Numbers indicate bootstrap support. CoVs are colored according to 


(Figure legend continued on the bottom of the next page.) 
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Box 2. Molecular Evolution in Viruses: Methods and Caveats 


Phylogenetic tree construction 

Distance-based and character-based approaches can be used in phylogenetic tree reconstruction. Distance-based 
methods measure pairwise differences among sequences and generate the tree from the resultant distance matrix (e.g., 
UPGMA, Neighbour Joining). Character-based methods evaluate all possible trees and estimate the one that best fits the 
data. This approach includes maximum likelihood (ML) and Bayesian methods (e.g., phyML [67], RaxML [68], MrBayes 
[69)). 


Recombination detection 
The location of recombination breakpoints can be detected using phylogenetic incongruence among segments in a 
sequence alignment (e.g., GARD [70}) or by evaluating the distribution of nucleotide substitution along genomic regions 
(e.g., Recco [71]). Acommon approach is to use different methodologies to identify breakpoint locations and contributing 
sequences (e.g., RDP [72)). 


Positive selection analyses 

Positive selection is usually estimated based on w variation across sites and/or lineages (Box 1). A common approach is 
to compare ML models that allow or not a class of codons in the alignment to evolve with w > 1 (e.g., the ‘site models’ in 
PAML [73]) Likelihood ratio tests are then applied to determine whether the neutral model can be rejected in favor of the 
positive selection model. Alternatively, branch-site models can be applied to detect episodic positive selection on specific 
branches of the phylogeny. Different methods allow to test a priori whether a branch is under selection [74,75] or to model 
different evolutionary scenarios for each branch [76]. When evolutionary analysis is focused on recent timescales, 
selection may not have yet fixed the advantageous mutations or removed the deleterious ones. One possibility is to 
compare the distribution of nonsynonymous and synonymous polymorphisms in a specific lineage with the ratio of 
nonsynonymous to synonymous fixed differences between lineages/species. The McDonald-Kreitman test [77] has been 
widely applied for this purpose. 


Relaxed selection 

When selection is relaxed, smaller w values tend toward 1, whereas w values higher than 1 decrease. This phenomenon 
can be confused with episodic selection. Specific methods allow one to infer whether a branch in the phylogeny is under 
positive or relaxed selection [18]. 


Subsitution rate saturation 

Saturation of substitution rates can be a serious issue for deep tree branches. Nonetheless, branch-site methods are 
relatively insensitive to biases introduced by dS saturation, and can be applied to the analysis of distantly related species 
[78]. Alternatively, specific indexes have been developed to detect substitution saturation [79]; in the presence of dS 
saturation, third-codon positions can be removed to obtain reliable phylogenies. 


Synonymous constraint 
Regions with an excess of synonymous constraint can be identified using a recently developed sliding-window ML-based 
method [80]. 


[14,15]. Analysis of the ORF8 region revealed high sequence identity with civet/numan SARS- 
CoV. Two groups came to the conclusion that recombination within SARSr-Rs-CoVs or 
between SARSr-Rs-CoVs and SARSr-Rf-CoVs led to the acquisition of an ORF8 closely related 
to that of civet/numan SARS-CoV and ultimately originated the virus responsible for the human 
epidemic [14,15]. Interestingly, Lau and coworkers also reported that the ORF8 region has a dN/ 
dS = 3.5 in SARS-CoVs isolated from humans (but not in SARSr-BatCoVs), indicating the action 
of positive selection (Box 1) [14]. This finding is interesting per se and becomes even more 
important considering that, early in the human epidemic, SARS-CoVs acquired a signature 29- 
nucleotide deletion which split ORF8 into two functional ORFs (ORF8a and b) [16]. These 
findings suggest that rapid evolution of ORF8 might facilitate host shifts [14]. This possibility is, 
however, questioned by the presence of additional SARS-CoV human isolates that carry 
independent and larger deletions in the ORF8 region [16]. Thus, an alternative explanation 


genus and lineage. Information about origin, intermediate host, and clinical presentation is reported for the six HCoVs [1-5,89]. Data about case fatality rate were derived 
from the World Health Organization website (htto://www.who.int/mediacentre/factsheets/mers-cov/; http://www.who.int/csr/sars/country/table2004_04_21/en/). (B) 
CoV genome organization is schematically reported together with information on receptor/coreceptor usage. Virus names are colored according to their genus or lineage, 
as in (A). Only ORFs mentioned in the text are colored or shaded. Empty boxes represent accessory ORFs that are not described in the text. 
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Box 3. Time Origin of CoV Genera and HCoV Emergence 


Coronaviruses are classified into four distinct genera (alpha, beta, gamma, and delta) [81]; aloha-CoVs and beta-CoVs circulate in mammalian hosts, whereas gamma- 
CoVs and delta-CoVs mainly infect birds [82]. An analysis of the RNA-dependent RNA polymerase (RdRp) gene in 43 CoVs provided a first estimate of around 10 000 
years ago for the time of the most recent common ancestor (tWRCA) of the four genera [83]. This result was questioned on the basis of observations suggesting a 
longstanding interaction between CoVs and their hosts [82]. Thus, Wertheim and coworkers hypothesized that natural selection, in particular negative selection, 
resulted in a bias of the tMRCA estimate [82]. Indeed, strong negative selection can result in the saturation of substitutions at synonymous sites and consequently in 
underestimation of branch lengths in a phylogeny. To overcome this issue, the authors applied a branch-site test (see Box 1) [76] to estimate branch lengths while 
taking into account the effect of different selective pressures among lineages in the CoV phylogeny [82]. Their findings placed the separation of the four CoV genera 
around 300 million years ago, highlighting the importance of evolutionary models in molecular clock dating. Interestingly, this tMRCA is consistent with the separation 
time between mammals and aves [84], suggesting a coevolutionary relationship between coronavirus and their hosts. However, the dating estimates obtained for 
closely related viruses by Wertheim and coworkers were in agreement with previous studies, suggesting that the action of natural selection is not biasing the estimation 
of more recent divergence times. Although most HCoVs were identified only recently, molecular clock analyses indicate that some of these viruses diverged from 
closely related CoVs hundreds of years ago. In particular, the emergence of HCoV-NL63 and HCoV-229E has been roughly estimated around 500-800 and 200 years 
ago, respectively [85,86]. HCoV-OC43 is thought to have shared a common ancestor with BCoV around 120 years ago [87]. As for SARS-CoV and MERS-CoV, 
molecular dating studies estimated that they diverged from bat CoVs in the last three decades [14,88]. Finally, the MRCA of HCoV-HKU1 extant lineages was estimated 
to have existed in the 1950s [60]. Clearly, these dates should be regarded as estimates and confidence intervals are often wide. A timeline for the emergence of HCoVs 
is depicted in Figure |. 
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Figure |. Timeline for the Emergence of HCoVs. 


for these findings is that the activity of ORF8 became dispensable in the human host. If this were 
the case, relaxed purifying selection rather than positive selection might be responsible for the 
high dN/dS. To disentangle these alternative possibilities we analyzed ORF8 in human and civet 
viruses that carry an intact gene, as well as in bat viruses. Although we confirmed that dN/dS is 
well above 1 for human/civet SARS-CoV ORF8, we detected no evidence of positive selection 
using the M7/MB ‘site models’ from PAML (Box 2) or with PARRIS (PARtitioning approach for 
Robust Inference of Selection) [17] (Figure 2A). Instead, we obtained evidence that relaxation of 
natural selection [18] in ORF8 accompanied the shift from bats to civets/numans (Figure 2A). 
These results suggest no major adaptive role for ORF8 during the human SARS-CoV epidemic 
and support the view that ORF8 is dispensable for virulence and transmission at least in the 
human/civet host. 


A similar gene loss from the genome of a bat-derived ancestor occurred during the evolution of 
HCoV-229E. CoVs closely related to HCoV-229E were recently isolated from African hippo- 
siderid bats [19], and a CoV belonging to the same species as HCoV-229E had been described 
in captive alpacas suffering from an acute respiratory syndrome [20,21] (Figure 1A). Analysis of 
these viral genomes indicated that, compared to HCoV-229E, they carry an additional ORF at 
the genomic 3’ end [20] (Figure 1B). This ORF, which is designated ORF8 but shares no 


Trends in Microbiology, Month Year, Vol. xx, No. yy 5 


TIMI 1370 No. of Pages 14 


Trends in Microbiology 


CelPress 


(A) 


Relaxed 2 
° 8 
oon. 2 
01 % 5 
os he 2 
1K Ss Zz . 
7 % Zz P value for relaxation (K=0) < 0.01 
z ro} 
% 2 
5 % . 
@ 
“aks o 
Intesified 


Civet SARS-CoV SZ3/2003 


Civet SARS-CoV SZ16/2003 


Sans Evolutionary analysis of complete ORF8 sequences 
e179! “Coy, 
espn cl S205 rack ~~ 1% 
ive! 
Number of sequences 6 18 
dN/dS (95% Cl) 0.18 (0.13, 0.25) 2.65 (1.33, 4.65) 
Test of positive selection: 
PAML site models Not performed P=0.18 
£ PARRIS Not performed P=0.67 
ro 
(B) Rat AKAP7 130 
MHV NS2a 46 
HCoV-0C43 NS2a 47 
Rotavirus A VP3 718 -4HLTLGVLFISDNM----IHDVITV...QSRMFI- 
MERS-CoV NS4b YH 
MERS-CoV NS4b MHV NS2a Rotavirus A VP3 
(Cc) 
SARS-CoV ORF3A 20 ------. IDNASPASTVHATA®?IPLQASLPFGWLV 1 GVArLAVFOSATKI IALNKRWOLALYKGFOF 1 CNL LLLP VTIYSHLLLVARGMEAQFL YLYAL! YFLOCINACRI IMRCWLCIKCKSKNPLLYDANYPV 147 
HCoV-229E ORF4a 27 EVSROVIOQDVKDGIVIFNLLAYTLMS-------- LFVVIFALFR)-----; -ARSHRGRAAL----- IVFRILILFV--YVPLLYWSQAYCYATLIAVILLG-------- REPHTAVHCWINKTWDFIVFNVTTLC 130 
HCoV-NL63 ORF3 28 HDVTVLRDNLKPVITLSTITAYLLVS-------- LFVIVFALFR)-----! PLTARGRVAC----— EVLELLTLFV--YVPLLVLFGMYLDSFIIFSTLLF-------- RETHVGIYAYLYKNFSFVLPNVTKLC 131 
TM1 TM2 TM3 
N-term N-term 
Extracellular 
Intracellular 
SARS-CoV ORF3A en 
HCoV-229E ORF4a 
HCoV-NL63 ORF3 HCoV-OC43 ORFS 
Trends in Microbiology 


(See figure legend on the bottom of the next page.) 
6 Trends in Microbiology, Month Year, Vol. xx, No. yy 


TIMI 1370 No. of Pages 14 


homology with the homonymous SARS-CoV gene, has unknown function and shows limited 
similarities to any other CoV gene [20]. We analyzed the sequences of recently identified alpha- 
CoVs from camels [22] and found that ORF8 is encoded by these viruses, as well (Figure 1B). 
Thus, it is presently unknown whether the loss of ORF8 conferred some advantage during the 
host shift to humans or, as in the case of ORF8 in SARS-CoV, it became dispensable in the 
human host. 


Another interesting feature of some CoVs is that they encode phosphodiesterases (PDEs) 
(Figure 1B). These viral enzymes cleave 2’,5’-oligoadenylate, the product of OAS proteins, to 
prevent activation of the cellular endoribonuclease RNase L and consequently block interferon 
(IFN)-induced antiviral responses [23]. The PDE activity in the mouse hepatitis virus (MHV) NS2a 
protein is critical for hepatovirulence [23]. HCoV-OC48, as well as other lineage A nonhuman 
beta-CoVs, encode NS2a proteins that are characterized by a high degree of sequence similarity 
to the MHV PDE (Figure 1B). A protein with structure and sequence homology to NS2a is also 
encoded by an unrelated virus, Group A rotavirus. In this case the PDE activity resides in the C- 
terminal portion of VP8, a virulence factor [24]. Interestingly, both VP3 and NS2a show two 
motifs that are characteristic of the 2H-PDE family and share very little sequence similarity to the 
PDE domain of a cellular protein, AKAP7 [24] (Figure 2B). AKAP7 and the viral PDEs display 
structural homology (Figure 2B), and murine AKAP7 can complement an inactive MHV NS2a 
gene [25]. From an evolutionary standpoint, these observations suggest that: (i) beta-CoVs and 
rotaviruses have independently acquired PDE activities; and (ii) AKAP7 served as the source 
gene in both viral genera (Figure 2B). More recently, a PDE activity was also discovered in the 
NS4b protein of MERS-CoV (Figure 1B) and other lineage C beta-CoVs [26]. Similar to those in 
lineage A beta-CoVs and rotavirus, NS4b belongs to the 2H-phosphoesterase family and 
displays a predicted structure homologous to AKAP7 [26] (Figure 2B). It remains to be deter- 
mined whether NS4b was acquired by capturing a vertebrate AKAP7, but the observation that 
distinct viruses acquired, most likely independently, a PDE activity underscores the importance 
of these enzymes for viral fitness. 


It was recently proposed that CoVs (and other viruses) stole additional genes from their hosts 
[27]. Hemagglutinin-esterases (HEs) are encoded by lineage A beta-CoVs (e.g., HCoV-HKU1 
and HCoV-OC48) (Figure 1B), as well as influenza C virus and toroviruses. Structural analysis 
suggested that these viral enzymes derive from an ancestral host lectin, although it is unclear 
whether acquisition occurred in an ancestral virus followed by speciation or multiple times [27]. 
Incidentally, the N-terminal domain of the CoV spike protein is also believed to derive from a 
cellular lectin [28]. Unlike the influenza virus C enzyme, CoV HEs lack membrane-fusion activity 
and are accessory to the spike protein by serving primarily as receptor-destroying enzymes 
(RDE) - that is, they aid viral detachment from carbohydrates present on infected cells [29,30]. In 
fact, HEs are present only in the genome of lineage A beta-CoVs, most of which use sialic acids 
as coreceptors [1] (Figure 1B). These observations suggest that sialic acid-binding spike 


Figure 2. Evolution of Human Coronavirus (HCoV) Accessory Proteins. (A) Test for relaxation of selective strength for SARS-CoV and SARSr-BatCoVs ORF8. 
Branches are colored according to the selection intensity parameter k. RELAX evaluates if selection on the test branches (bold) is relaxed (k < 1) or intensified (k > 1) 
compared to background branches. In the evolutionary analysis table the number of sequences differs from that in the tree because RELAX removes identical sequences. 
Evidence of positive selection was searched for using the M7/M8 ‘site models’ from PAML (see Box 2) and with PARRIS. M7 and M68 represent the null and the positive 
selection models, respectively. A likelihood ratio test (with 2 degrees of freedom) was applied. (B) An amino acid alignment of rodent AKAP7 and four viral 
phosphodiesterases (PDEs) is shown. Amino acids are colored red if they are identical, orange if they have very similar properties. PDEs belonging to the 2H family 
are characterized by two H-®-[S/T]-® motifs (blue boxes), where ® is a hydrophobic residue. The structure of rat AKAP7 (gray, PDB ID: 2VFK) is superimposed on MERS- 
CoV NS4b (green, model generated from 2VFK), MHV NS2a (cyan, PDB ID: 4Z5V), and Rotavirus A VPS (yellow, PDB ID: SAF2). Catalytic histidines are shown in red. (C) 
Sequence and membrane topology comparison of HCoV viroporins. Transmembrane regions (TM1-3) predicted by the TMHMM algorithm [90] are boxed in blue. The 
corresponding topology model for SARS-CoV ORF3A, HCoV-229E ORF4a (from the Inf-1 strain), and HCoV-NL63 ORFS is shown. The topology model of HCoV-OC43 
OFR5 was derived from recent data [34]. 
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proteins coevolved with HE genes serving as RDEs. This hypothesis is supported by the 
observation that the MHV spike protein evolved from an ancestral sugar-binding preference 
to a protein-binding mode and that several MHV strains lost expression of HE [27,28] (Figure 1B). 


Finally, it is important to notice that artificial selection can lead to unintended changes in viral 
genomes. Such changes most likely result from passages in culture that, on one hand, relieve the 
virus from pressures exerted in vivo (e.g., by the host immune system) and, on the other hand, 
derive from viral adaptation to the in vitro system. An example of these effects is the loss of a full- 
length ORF4 in the HCoV-229E prototype strain which, due to a two-nucleotide deletion, has a 
split gene, encoding two proteins (ORF4a and ORF4b) [31,32] (Figure 1B). Conversely, clinical 
isolates display a full-length ORF4 sequence [32]. An intact ORF4 is also observed in bat and 
camel viruses related to HCoV-229E [19,22], whereas the alpaca alpha-CoV displays a one- 
nucleotide insertion, resulting in a frameshift [20] (Figure 1B). The availability of only a single 
alpaca CoV genome makes it impossible to determine whether the inserted sequence is 
representative of the alpaca CoV population or, else, if it represents a sequencing error. 


Overall, these observations suggest that loss of full-length ORF4 is a result of passaging in cell 
culture, a process that often generates attenuated viruses. An interesting finding on the role of 
ORF4a came from the observation that its protein product regulates virus production in vitro by 
functioning as a viroporin [33]. Most likely, the same function is performed by the full-length 
ORF4 as well. 


Viroporins were also detected in SARS-CoV, HCoV-OC43, and HCoV-NL68 [34,35] (Figure 1B). 
As expected, given the relatedness of the two viruses, the proteins from HCoV-NL63 and HCoV- 
229E share substantial sequence similarity. Limited similarity is also observed with the SARS- 
CoV protein, especially in the first and second transmembrane regions, suggesting either a 
common origin or independent acquisition followed by convergent optimization of residues in the 
transmembrane helices (Figure 2C). Conversely, the HCoV-OC43 protein (encoded by ORF5, 
originally denoted NS12.9) is unrelated to the other CoV viroporins, both in terms of sequence 
and of domain topology [34] (Figure 2C). A protein homologous to the HCoV-OC48 viroporin is 
instead encoded by MHV (accessory protein NS5a) and functions as an antagonist of IFN- 
induced antiviral responses [34,36]. Whether the HCoV-OC48 viroporin has the same IFN- 
antagonizing activity remains to be investigated; however, mutant viruses lacking ORF5 display 
growth defects in vitro and in vivo, as well as reduced virulence in mice [34]. Interestingly, the 
viroporins from SARS-CoV, HCoV-NL63, and HCoV-229E can complement the viroporin- 
defective mutant HCoV-OC48 in vitro [34]. Thus, the conserved function of CoV viroporins 
was most likely attained by convergent evolution for acquisition of unrelated genes. 


Evolution of Structural and Nonstructural Proteins 
Clearly, CoV genomes do not only evolve by gene gains and losses, but also via subtler changes 
that modify protein sequences, and recombination has an important role in reassorting variants. 


£ 


Several excellent reviews have focused on the evolutionary history of SARS-CoV genomes in 
terms of recombination and natural selection [37-39]; hereafter, SARS-CoV will be mentioned 
only to draw comparisons with other CoVs. 


From an evolutionary standpoint, nonstructural proteins have attracted less attention than the 
structural components. This is likely due to the fact that proteins exposed on the virus surface 
represent the preferential targets of the host immune response. Thus, analyzing and describing 
their variability and evolutionary dynamics has a clear relevance for the development of preven- 
tive strategies (e.g., vaccines) and of treatment options (e.g., administration of neutralizing 
antibodies). Moreover, structural proteins, and the S protein in particular, determine the first 
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and essential steps in infection and most likely represent the major determinants of host and 
tissue tropism. 


In CoVs, the S protein includes two functionally distinct units: the S1 region contains an N- 
terminal domain (NTD) and the receptor-binding domain (RBD, also referred to as C-terminal 
domain or CTD), whereas the S2 region includes the fusion peptide, two heptad repeats (HR1 
and HR2), and the transmembrane region (Figure 3A) [38]. A striking feature of HCoV spike 
proteins is that they have adapted to use diverse cellular receptors and there is no congruence in 
the phylogeny of HCoV and their receptor usage. In fact, closely related viruses may use different 
receptors (Figure 1B). For instance, HCoV-229E uses aminopeptidase N (ANPEP), whereas 
HCoV-NL63 exploits ACE2, this latter representing the receptor for the relatively divergent 
SARS-CoV (Figure 1B). It is presently unclear how these binding specificities evolved. The latest 
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Figure 8. Evolution at the Coronavirus (CoV)—Host Interaction Surface. (A) Schematic representation of MERS-CoV 
spike protein domains. Positively selected sites in MERS-CoV and other lineage C beta-CoVs are shown in red, RBD 
mutations emerged in the South Korean outbreak are in magenta (see text). A detail of the interaction surface between the 
MERS-CoV RBD and human DPP4 (PDB ID: 4F5C) is also reported. (B) Ribbon diagram of the interaction surface of human 
ACE2 with the spike protein of SARS-CoV (PDB ID: 2AJF) and HCoV-NL63 (PDB ID: 83KBH). The binding surface of porcine 
ANPEP with the TGEV spike protein (PDB ID: 4F5C) is also shown. The location of the HCoV-299E binding site on ANPEP is 
circled. Red denotes protein regions involved in binding. 
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developments on this topic and, more generally, on the evolution of structural and nonstructural 
proteins are detailed below for the five known HCoVs. 


MERS-CoV 

The evolutionary analysis of MERS-CoV is a rapidly moving field, as sequences from the latest 
phases of the epidemic have just become available. Analysis of an ever increasing number of viral 
sequences of both MERS-CoV and of related beta-CoVs revealed that genetic variability in the S 
gene was shaped by recombination and positive selection. In fact, both ancient and recent intra- 
spike recombination events were described [22,40,41]. Interestingly, recombination events with 
breakpoints within the S gene occurred in camels in Saudi Arabia and originated the MERS-CoV 
lineage that spread to South Korea. 


Analysis of positive selection of MERS-CoV spike genes indicated that several adaptive variants 
arose in MERS-CoV and in phylogenetically related CoVs [42]. Contrary to common expectation 
and to what happened during the SARS-CoV host shift to humans, positive selection did not 
target the RBD. In fact, most adaptive substitutions were detected in the region encompassing 
the heptad repeats, regions of central importance for virus cell entry (Figure 3A) [42,43]. In other 
CoVs, variants in the heptad repeats were previously shown to affect host or tissue tropism [44— 
46]. Interestingly, during the South Korean outbreak, MERS-CoVs that carry point mutations in 
the spike protein RBD emerged and rapidly spread [47]. These viruses showed decreased 
binding to the cellular receptor [47] (Figure 3A). Because several immune epitopes are located in 
the RBD, these findings point to the possibility that MERS-CoV is evolving to avoid the binding of 
neutralizing antibodies, resulting in a trade-off with receptor-binding affinity [47]. If this were the 
case, the phases of MERS-CoV adaptation to humans may have consisted of initial events that 
modulated host tropism through changes in the heptad repeats followed by the emergence of 
virus variants that escape immune responses. In MERS-CoV and other lineage C beta-CoVs, 
positive selection also targeted the nonstructural components, particularly ORF1a [48]. Most 
adaptive events occurred in nsp3, a multifunctional protein which acts as a viral protease and 
contributes to the suppression of interferon responses through its deubiquitinating and del- 
SGylating activities [49]. Selection in nsp3 is ongoing among MERS-CoV isolated from humans 
and camels [48]. In analogy to the S protein, though, no major selective event was found to be 
associated with camel-to-human transmission, although a positively selected change (R911C) in 
nsp3 was observed among human-derived viruses alone, suggesting that viral adaptation to our 
species represented the underlying pressure [48]. 


HCoV-229E 

A recent analysis indicated that HCoV-229E may have recombined with the alpaca alpha-CoV 
virus within the S gene, as also demonstrated by the distinct phylogenetic trees for the S1 and S2 
regions [19]. Also, HCoV-229E acquired a deletion in the S gene compared to bat viruses [19]. 
Recent sequencing of several of such viruses showed that this deletion is also present in the 
alpaca CoV S gene and in camel-derived alpha-CoVs [22]. This finding is particularly interesting 
because deletions in the NTD are associated with changes in tissue tropism in TGEV (trans- 
missible gastroenteritis virus): in this porcine virus the spike has dual tropism for the respiratory 
and intestinal tracts, but the N-terminally deleted variants from PRCV (porcine respiratory 
coronavirus) only infect the respiratory tract [50,51]. In chiroptera, CoVs are mainly restricted 
to the gastrointestinal tract, whereas in humans and camelids, the upper and lower respiratory 
airways are infected. It will be important to determine whether the S gene deletion in HCoV-229E 
and camelid alphaCoVs is indeed responsible for a change in tissue tropism. 


HCoV-NL63 


Recombination contributed to shaping the diversity of the S gene among HCoV-NL63 viruses. 
Recombination between an ancestral HCoV-NL63 virus and the related PEDV was also 
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detected in the M gene that, in its 3’ portion, is more similar to PEDV than to HCoV-229E [52]. 
Like SARS-CoV, HCoV-NL63 uses its RBD to bind ACE2. The binding site on the cellular 
receptor is the same for the two viruses but the RDBs show no sequence similarity. Interestingly, 
the RBDs of SARS-CoV and HCoV-NL68 do not display any structural similarity either: HCoV- 
NL63 contacts ACE2 with three discontinuous beta-loops, whereas SARS-CoV binds the 
receptor through a continuous subdomain [53] (Figure 3B). These observations suggest that 
the two viruses independently acquired the ability to bind the same ACE2 region via convergent 
evolution or that they shared an ACE2-binding ancestor long ago. Strikingly, TGEV, which is 
phylogenetically related to HCoV-NL63, uses two regions corresponding to the HCoV-NL63 
beta-loops to bind a distinct cellular receptor, ANPEP (Figure 1B, Figure 3B) [54]. Finally, HCoV- 
229E, sharing sequence homology with HCoV-NL63 and TGEV (Figure 1A), binds ANPEP, but 
engages a region distinct from that bound by TGEV [55]. Overall, these data highlight the 
extraordinary plasticity of CoV RBDs, and their complex evolutionary dynamics whereby diver- 
gent evolution can be followed by convergent adaptation to the same receptor. This complexity 
is further expanded by the ability of some CoVs to use other cellular attachment molecules to 
complement the function of the RBD. Indeed, the S protein of HCoV-NL63 exploits heparan 
sulfate proteoglycans to adhere to host cells [56]. Interestingly, a similar ability to bind heparan 
sulfate can be gained by MHV with relatively few in vitro-acquired mutations in the S protein [57]. 
In line with the view that heparan sulfate is an aspecific receptor, the mutant MHV viruses display 
expanded host tropism [57], highlighting the potential relevance of combinatorial receptor usage 
or receptor shifts for interspecies transmission. 


HCoV-OC43 

Recombination seems to be rampant in HCoV-OC48 viruses and contributed to originate the A 
to E viral genotypes, as well as viruses that do not belong to these major genotypes [58-60]. To 
our Knowledge, no study has analyzed the fitness of recombinant viruses or, more generally, of 
viruses belonging to distinct genotypes. Nonetheless, two reports indicated that genotype D has 
become predominant in the East Asian population [58,59]. Whether this is due to population 
acquired immunity against the older A and B genotypes or to viral features unrelated to 
antigenicity remains to be determined. 


The active recombination in HCoV-OC43 suggests that inference of natural selection is best 
performed by analysis of sequences belonging to the same genotype. In one such analysis, 
positive selection was found to act on the S gene of genotype D viruses [61]. Interestingly, 
several positively selected sites with high posterior probability of positive selection are located in 
the NTD, which is involved in the binding of sialic acids. 


A positively selected site was located in the CTD, a region that has unclear function in the HCoV- 
OC48 S protein, as no known protein receptor has been identified to date [61]. However, recent 
data from HKU1 suggest that, by analogy, a protein receptor for HCoV-OC43 may exist [62] (see 
below). 


HCoV-HKU1 

The structure of the S protein of HKU1 was recently solved; the glycan-binding site is located in 
the NTD and is conserved with bovine coronavirus (BCoV) S1 [63]. Nonetheless, antibodies 
against the CTD, but not those against the NTD, block HKU1 infection of human tracheal- 
bronchial epithelial cells, suggesting that the CTD is the major RBD, and that a protein receptor 
for HKU1 exists [62]. In analogy to HCoV-NL63, glycans may mediate only the initial attachment 
to the host cells. 


A recent survey of HKU1 clinical isolates from different geographic origins indicated that most 
viruses from Colorado form a subclade in the HKU1 phylogeny and carry three distinctive 
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substitutions in the S protein within the NTD, CTD, and close to the S1/S2 cleavage site (W197F, 
F613Y, and H716D, respectively) [64]. It will be interesting to assess whether these differences 
are functional and derive from a selective process. 


Concluding Remarks 

Thanks to high-throughput techniques, a large number of complete CoV genomes have become 
available to the scientific community, and many more will be coming in the near future. Field 
studies have contributed enormously to widen our knowledge on the diversity of CoVs hosted by 
different vertebrates, and epidemiological surveys have provided CoV sequences from distinct 
geographic areas and associated with different disease phenotypes. In parallel, resources have 
been created to store and mine these data (e.g., The Virus Pathogen Database and Analysis 
Resource, VIPR [65]). These advances have allowed tracing the evolutionary history of the large 
and complex CoV genomes to an unprecedented detail. The emerging picture indicates that 
CoV genomes display high plasticity in terms of gene content and recombination. The long CoV 
genome expands the sequence space available for adaptive mutation, and the spike protein can 
adapt with relative ease to exploit different cellular receptors. These features are likely to underlie 
the alarming propensity of CoVs for host jumps. Despite these advances, major challenges 
remain (see Outstanding Questions). Efforts to underscore the viral genetic determinants that 
favor interspecies transmission should be pursued as an effective strategy to prevent or prepare 
for future HCoV emergence. 


Supplemental Information 
Supplemental information associated with this article can be found online at doi:10.1016/j.tim.2016.09.001. 


References 

1. Graham, R.L. et al, (2013) A decade after SARS: strategies for 3. Yang, X.L. et al. (2015) Isolation and characterization of a novel bat 
controlling emerging coronaviruses. Nat. Rev. Microbiol. 11, 836- coronavirus closely related to the direct progenitor of severe acute 
848 respiratory syndrome coronavirus. J. Virol. 90, 3253-3256 

2. Lau, S.K. et al. (2015) Discovery of a novel coronavirus, China 4. Lau, S.K. et al. (2015) Severe acute respiratory syndrome (SARS) 
Rattus coronavirus HKU24, from Norway rats supports the murine coronavirus ORF8 protein is acquired from sars-related coronavi- 
origin of Betacoronavirus 1 and has implications for the ancestor of rus from greater horseshoe bats through recombination. J. Virol. 
Betacoronavirus lineage A. J. Virol. 89, 3076-3092 89, 10532-10547 

3. Chan, J.F. et al. (2013) Interspecies transmission and emergence 5. Wu, Z. et al. (2016) ORF8-related genetic evidence for Chinese 
of novel viruses: lessons from bats and birds. Trends Microbiol. horseshoe bats as the source of human severe acute respiratory 
21, 544-555 syndrome coronavirus. J. Infect. Dis. 213, 579-583 

4. Drexler, J.F. et al. (2014) Ecology, evolution and classification of 6. Chinese SARS Molecular Epidemiology Consortium (2004) Molec- 


bat coronaviruses in the aftermath of SARS. Antiviral Res. 101, ular evolution of the SARS coronavirus during the course of the 
45-56 SARS epidemic in China. Science 303, 1666-1669 

5. Su, S. et al. (2016) Epidemiology, genetic recombination, and 7. Scheffler, K. et al. (2006) Robust inference of positive selection 
pathogenesis of coronaviruses. Trends Microbiol. 24, 490-502 from recombining coding sequences. Bioinformatics 22, 2493- 

6. Eckerle, L.D. et al, (2007) High fidelity of murine hepatitis virus 2499 
replication is decreased in nsp14 exoribonuclease mutants. J. 8. Wertheim, J.O. et al. (2015) RELAX: detecting relaxed selection in 
Virol. 81, 12135-12144 a phylogenetic framework. Mol. Biol. Evol. 32, 820-832 

7. Eckerle, L.D. et al. (2010) Infidelity of SARS-CoV Nsp14-exonu- 9. Corman, V.M. et al. (2015) Evidence for an ancestral association of 
clease mutant virus replication is revealed by complete genome human coronavirus 229E with bats. J. Virol. 89, 11858-11870 


sequencing. PLoS Pathog. 6, e1000896 20. Crossley, B.M. et al. (2012) Identification and characterization of a 
8. Vega, V.B. et al. (2004) Mutational dynamics of the SARS coro- novel alpaca respiratory coronavirus most closely related to the 
navirus in cell culture and human populations isolated in 2003. human coronavirus 229E. Viruses 4, 3689-3700 


BMC Infect. Dis. 4, 32 21. Crossley, B.M. et al. (2010) Identification of a novel coronavirus 

9. Lauber, C. et al. (2013) The footprint of genome architecture in the possibly associated with acute respiratory syndrome in alpacas 
largest genome expansion in RNA viruses. PLoS Pathog. 9, (Vicugna pacos) in California, 2007. J. Vet. Diagn. Invest. 22, 94— 
e1003500 97 


10. Subissi, L. et al. (2014) One severe acute respiratory syndrome 22. Sabir, J.S. et al. (2016) Co-circulation of three camel coronavirus 
coronavirus protein complex integrates processive RNA polymer- species and recombination of MERS-CoVs in Saudi Arabia. Sci- 
ase and exonuclease activities. Proc. Natl. Acad. Sci. U. S.A. 111, ence 351, 81-84 
E3900-E3909 23. Zhao, L. et al. (2012) Antagonism of the interferon-induced OAS- 

11. Menachery, V.D. et a/. (2016) SARS-like WIV1-CoV poised for RNase L pathway by murine coronavirus ns2 protein is required for 
human emergence. Proc. Natl. Acad. Sci. U. S.A. 113, 3048- virus replication and liver pathology. Cell. Host Microbe 11, 607- 
3053 616 

12. Ge, X.Y. et al. (2013) Isolation and characterization of abat SARS- 24. Zhang, R. et al. (2013) Homologous 2’,5’-phosphodiesterases 
like coronavirus that uses the ACE2 receptor. Nature 503, 535- from disparate RNA viruses antagonize antiviral innate immunity. 
538 Proc. Natl. Acad. Sci. U. S. A. 110, 18114-13119 


12 Trends in Microbiology, Month Year, Vol. xx, No. yy 


CelPress 


Outstanding Questions 


Do different genotypes (arising from 
recombination or mutation) of the same 
CoV species result in distinct viral 
phenotypes? 


What is the relevance of adaptive 
changes in nonstructural proteins for 
host adaptation? And what is the role 
played by accessory proteins in inter- 
species transmission and virulence? 


Why do similar viruses determine very 
different disease phenotypes in distinct 
mammalian hosts? Do the viral and 
host genome interplay in determining 
disease? And to what extent? 


What is the distribution of CoVs in dif- 
ferent mammalian orders? And how 
genetically diverse are CoVs hosted 
in wild or domestic mammals? 


Can evolutionary information be 
applied to develop tools to predict 
which viruses have greater zoonotic 
potential and may result in the next 
epidemic? 
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